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f *^ ' Abstract. A Noetherianity criterion is provided in this paper for the quo- 

r*f\ , tients of the algebra of partial difference polynomials with respect to its differ- 

ence ideals. After we introduce a Grobner bases theory for ideals generated by 
difference polynomials, this implies the existence of finite such bases for the dif- 
rj ' ference ideals containing elements with suitable linear leading monomials. As 

^Ji, , a by-product, one obtains computational methods for ideals of finitely gener- 

^~> i ated polynomial rings that are invariant under the action of finite dimensional 

pH ■ commutative algebras and in particular of finite abelian groups. 

-(— > ' 

1. Introduction 

The theory of difference algebras (see books [7] [H] [22] and references therein) 
was introduced in the 1930s by the mathematician Joseph Fels Ritt at the same 
time of the theory of differential algebras. Indeed, for a quite long time, difference 

C^ ' algebras has attracted less interest among researchers in comparison with differen- 

tial ones despite the fact that numerical integration of differential equations relies 
on solving finite difference equations. The rapid development of symbolic compu- 
tation and computer algebra in the last decade of the previous century gave rise to 
rather intensive algorithmic research in differential algebras and to creation of so- 
phisticated software as the diffalg library [3], implementing the Rosenfeld- Grobner 
algorithm and included in Maple and the package LDA [13]. At the same time, 
except for algorithmization and implementation in Maple of the shift algebra of 
linear operators [6] as a part of the package Ore_ALGEBRA, practically nothing has 

^\. , been developed in computer algebra in relation to difference algebras. Nevertheless, 

in the last few years, the number of applications of the theory and the methods of 
difference algebras is increased fastly. For instance, it turned out that difference 
Grobner bases may provide a very useful algorithmic tool for reduction of multi- 
loop Feynman integrals in high energy physics 9 , for automatic generation of finite 
difference approximations to partial differential equations [TTJ [5T] and for the con- 
sistency analysis of these approximations |1Q[ \T2\ . Another source of applications is 
related to the notion of "letterplace correspondence" [16l [17] [19] which transforms 
non-commutative computations for presented groups and algebras into analogue 
computations for ordinary difference algebras. As a result of all this use, a number 
of computer algebra packages implementing involutive and Buchberger's algorithms 
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2 V. GERDT AND R. LA SCALA 

for computing difference Grobner bases has been developed (see [101 Q2J [18] and 
reference therein). A major drawback in these computations, as for the differential 
case, is that such bases may be infinite owing to Non-Noetherianity of the algebra 
of difference polynomials. In fact, if X is a finite set and E denotes a multiplicative 
monoid isomorphic to (N r , +) then the algebra of difference polynomials is by defi- 
nition the polynomial ring P in the infinite set of variables IxS. Then, to provide 
the termination of the algorithms computing Grobner bases in P at least in some 
significant cases, we propose in this paper essentially two solutions. One consists 
in defining an appropriate grading for P that allows finite truncated computations 
for difference ideals J C P generated by a finite number of homogeneous elements. 
For monomial orderings of P that are compatible with such grading this implies 
a criterion, valid also for the non-graded case, able to certify the completeness of 
a finite Grobner basis computed on a finite number of variables of P. After the 
algebra of partial difference polynomials and its Grobner bases are introduced in 
Section 2 and 3, this approach is described in Section 4 and an illustrative exam- 
ple based on the approximation of the Navier-Stokes equations is given in Section 
5. A second solution to the termination problem consists in requiring that the 
difference ideal J contains elements with suitable linear leading monomials which 
corresponds to have the Noetherian property for the quotient algebra P/J. Some 
similar ideas appeared for the differential case in [5J [21] • One finds this second 
approach in Section 6. It is interesting to note that an important class of such 
Noetherian quotient algebras is given by polynomial rings P' in a finite number 
of variables which are under the action of a tensor product of a finite number of 
finite dimensional algebras generated by single elements. These finite dimensional 
commutative algebras includes, for instance, group algebras of finite abelian groups 
and hence, as a by-product of the theory of difference Grobner bases, one obtains 
a theory for Grobner bases of ideals of P' that are invariant under the action of 
such groups or algebras (see also [HI El]). These ideas are presented in Section 7 
and a simple application is described in Section 8. Finally, in Section 9 one finds 
conclusions and hints for further developments of this research. 



2. Algebras of difference polynomials 

Let (E, •) be a monoid with identity element I and let if be a field. We denote 
by End(if ) the monoid of ring endomomorphisms of K . We say that E acts on 
K or equivalently that if is a E- field if there exists a monoid homomorphism 
p : E -> End(if). In this case, we denote a ■ a — p(a)(a), for all a G E and 
a G K . Starting from now, we assume that if is a E-field. Let A be a commutative 
if-algebra. We say that A is a Yj-algebra if there is a monoid homomorphism 
p' : E -> End(yl) extending p : E -> End(if) that is p'{(j){a) = p{o-)(a) for all 
a G E and a G K . Let B C A be a if-subalgebra of a E-algebra A. We call 
B a E-subalgebra of A if E • B C B. In the same way, one defines E-ideals of 
A. Let B G A be a E-subalgebra and let X c B be a subset. We say that B is 
T,-generated by X if B is generated by E • X as a if-subalgebra. In other words, B 
coincides with the smallest E-subalgebra of A containing X. In the same way, one 
defines E-generation for E-ideals. Let A, B be E-algebras and let ip : A — > B be a 
if-algebra homomorphism. We call (p a E -homomorphism if f(a ■ x) = a ■ (p(x) for 
all a G E, x G A. 
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In the category of E-algebras one can define free objects as follows. Let X be a 
set and denote by x(a) each element (x, a) of the product set X(E) =XxE. Define 
P = K[X(E)] the polynomial algebra in all commuting variables x(cr) e -X"(S). For 
any element a £ E consider the E-algebra endomorphism a : P — > P such that 
x(t) i-> x(ar), for all x(t) £ X(Ei). Then, one has a faithful monoid representation 
p : E — > End(P) such that p(a) — a ad hence P is a E-algebra. Note that if E is a 
left-cancellative monoid then all maps p(cr) are injective. 

Proposition 2.1. Let A be a "E-algebra and let f : X — > A be any map. Then, a 
unique E-algebra homomorphism ip : P — > A is given such that tp(x(l)) = f{x), for 
all x £ X . 

Proof. A i^-algebra homomorphism ip : P — >• A is clearly defined by putting 
<p(x(o~)) = a ■ f(x), for any x £ X and <j S E. Then, one has that p(j ■ ax (a)) = 
p((t ■ a)x(ra)) — (r • a)p(x{ra)) — (r ■ a)(ra ■ f(x)) — r • (a(o~ ■ f(x))) — 
t ■ (ap(x(o~))) — t ■ tp(ax(o~)), for all a £ K, x G X and a, r G E. □ 

Definition 2.2. We call P — K[X(T,)] the free E-algebra generated by X . In fact, 
P is Yi-generated by the subset X(l) = {xi(l) \ Xi £ X}. 

Then, if A is any E-algebra which is E-generated by X one has that A is iso- 
morphic to the quotient P/J where J C P is the E-ideal containing all E-algebra 
relations satisfied by the elements of X . In other words, there is a surjective E- 
algebra homomorphism ip : P — > A such that x(l) H- x (x € X) and one defines 
J = Ker p. 

From now on, we specialize our assumptions in the following way. Let X = 
{xi, . . . , x n } be a finite set and let E be a free commutative monoid generated by 
a finite set, say {o~i, . . . ,oy}. Note that (E, •) is a cancellative monoid isomorphic 
to (N r ,+). Under the above assumptions, the finitely generated free E-algebra 
P = K[X(T,)] is called the algebra of partial difference polynomials. In literature 
one finds also the notation P = K{X} which emphasizes the role of X as (free) 
E-generating set of the algebra P. If if is a field of constants that is E acts 
trivially on K then the difference polynomials of P are with constant coefficients. 
Moreover, one uses the term ordinary difference when r = 1. The motivation 
for these names is the following. Let K be a field of functions in the variables 
t\, . . . , t r and fix hi, . . . , h r some parameters (mesh steps). Assume we may define 
the action of E on K by putting a ■ f(t\, . . . , t r ) = f(t\ +a\hx, . . . , t r + a r h r ) G K 
for all a — FJi a T 6 ^ an d for any function / £ K. For instance, one can 
consider the field of rational functions K — F(t\, . . .,tk) over some field F and 
hi, . . . , h r G F. Then, in the formal (symbolic) theory of difference equations the 
variables Xi(l) are by definition unknown functions itj(ii, . . . , t r ) that are considered 
algebraically independent together with the shifted functions Xi{o~) = a • Xj(l) = 
Uiiti +aihi, . . . ,t r + a r h r ). Then, a E-ideal / C P is also called a partial difference 
ideal and a E-basis of / corresponds to a system of partial difference equations. Since 
E is an infinite set, we have that P = K[X (E)] is not a Noetherian ring and hence 
such ideals have bases or E-bases which are generally infinite. 

3. Grobner bases of difference ideals 

In this section we introduce a Grobner basis theory for the algebra of partial 
difference polynomials by extending what has be done in [18] for the case of constant 
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coefficients. Note that the concept of difference Grobner basis has arisen also in 
[TO) IT3J [T7] . We start by noting that the Higman's Lemma [M] implies that the 
polynomial algebra P can be endowed with a monomial ordering even if the variables 
set -XT(E) is infinite. Denote by M = Mon(P) the set of all monomials of P. 

Proposition 3.1. Let ~< be a total ordering on M such that 

(i) 1 -< m, for all m G M, m =/= 1; 
(ii) if rn -<n then tm -< tn, for any m, n, t G M ; 
(iii) -< is a well- ordering on -X"(E) C M . 
Then, ■< is a well- ordering also on M that is -< is a monomial ordering on P. 

Note that the monoid E stabilizes the monomials set M. Then, it is clear that 
one needs monomial orderings that are compatible with the action of E on M. 

Definition 3.2. Let -<, be a monomial ordering of P. We call -< a monomial 
E-ordering of P if m ■< n implies that a ■ m -< a ■ n, for all m,n G M and a G E. 

Note that if -< is a monomial E-ordering of P then one has immediately that 
a ■ m >z rn, for all m G M and a G E. Examples of such orderings can be easily 
constructed in the following way. Let Q = K [a\, . . . , ay] be the polynomial algebra 
in the variables o~j and therefore E = Mon(Q). Moreover, let K[X] = K[xi, . . . , x n ] 
be the polynomial algebra in the variables Xi . Fix a monomial ordering < for Q and 
a monomial ordering -< for ii"[X]. For any a G E denote X(a) = {xi(a) x% G X}. 
Clearly P(cr) = K[X(a)] is a subalgebra of P that is isomorphic to K [X] and hence 
it can be endowed with the monomial ordering -<,. Since P = (^ ffeS P{°~) one can 
define a block monomial ordering for P obtained by < and -<. 

Proposition 3.3. Let m = m(5i) ■ ■ •m(<5/ 8 ),n = n(Si) ■ ■ -n(5k) G M where m(Si), 
n(5i) G M(Si) = Mon(P((5i)) and Si > . . . > Sk- We define m -<! n if and only if 
there is 1 < i < k such that m(Sj) = n(5j) when j < i and m(Si) -< n(8i). Then, 
-<' is a monomial Y,-ordering of P. 

Proof. For all a G E, one has that a ■ m = m(o~5i) ■ ■ ■ m(aSk) where m(o~5i) — 
a ■ m(6i) G M(o~8i) and a5\ > . . . > aSk since < is a monomial ordering of Q. 
Assume m <' n that is m(Sj) = n(Sj) for j < i and m(5i) -< n(Si). Clearly, one has 
also that m(o~Sj) — n(o~5j). Moreover, by definition of the monomial ordering -<, on 
all subalgebras P{o~) C P we have that m(5i) ~< n(8i) if and only if m(ab~i) -i, n(a8i). 
We conclude that a ■ m -<' a ■ n. D 

From now on, we assume that P is endowed with a monomial E-ordering -<. Let 
/ = J2i c i m i 6 P w ith m i e ^^ an d 7^ Ci G K. We denote as usual lm(/) = 
nik = max^{rrii}, lc(/) = c^ and lt(/) = lc(/)lm(/). Since -< is a E-ordering, one 
has that lm(cr ■ /) = a ■ lm(/) and therefore lc(<7 • /) = a ■ lc(/), lt(cr • /) = a ■ lt(/), 
for all a G E. If G G P we put lm(G) = {hn(/) / G G, f ^ 0} and we denote as 
LM(G) the ideal of P generated by lm(G). 

Proposition 3.4. Let G C P. Then lm(E • G) = E ■ lm(G). In particular, if I is 
a Yi-ideal of P then LM(/) is also H-ideal. 

Proof. Since P is endowed with a E-ordering, one has that lm(cr ■ /) = a ■ hn(/) for 
any /eP,//0 and a G E. Then, E • lm(I) = lm(E • I) C lm(I) and therefore 
LM(J) = (lm(/)} is a E-ideal. D 
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Definition 3.5. Let I C P be a H- ideal and G C I. We call G a Grobner E-basis 
of I iflm(G) is a Y>-basis o/LM(7). In other words, £ • G is a Grobner basis of I 
as an ideal of P. 

Since P is not a Noetherian algebra, in general its S-ideals have infinite (Grobner) 
£-bases. Note that one has a similar situation for the free associative algebra and 
its ideals and this case is strictly related with the algebra of ordinary difference 
polynomials owing to the notion of "letterplace correspondence" [HH [L7\ [19] . Nev- 
ertheless, in Section 6 we will prove the existence of a class of £-ideals containing 
finite Grobner E-bases. According to [101 E], such finite bases are also called 
"difference Grobner bases" . 

Let f,g £ P, f,g ^ and put lt(/) = cm, lt(g) = dn with m,n G M and c,d G K. 
If I = lcm(m,n) one defines the S-polynomial spoly(/, g) = (l/cm)f — (l/dn)g. 

Proposition 3.6. For all f.g £ P,f,g ^ and for any a £ £ one has that 
o- ■ spoly (/, g) = spoly(er • f,a ■ g). 

Proof. Note that lt(er-J) = (a ■ c) (a ■ m) , It (a ■ g) = (a-d)(o-n) with a-m,o-n £ M 
and a ■ c, o ■ d £ K . Since £ acts on the variables set X(S) by injective maps, if 
I = lcm(m, n) then a ■ I — lcm(i7 ■ m,a ■ n) and therefore we have 

a-spoly(f,g) = a-(^f-^g) = 
(*-c)(i- m ) a ■ f - (aI)(ln) (T ■ 9 = spoly(a • f,a ■ g). 

□ 
In the theory of Grobner bases one has the following important notion. 

Definition 3.7. Let f £ P, f ^ and G C P. If f = J2 t fi9i mth fi € p ,9i € G 
and lm(/) y lm(/j)lm((7j) for all i, we say that f has a Grobner representation 
with respect to G. 

Note that if / = ^- ftgi is a Grobner representation then a- f = J2i( a ' fi)( a '9i) 
is also a Grobner representation, for any o £ E. In fact, from lm(/) y ]xa(fi)lax(gi) 
it follows that lm(cr • /) = a ■ hn(/) y (a ■ lm(/j))((T • lm(pj)) = lm(cr • /i)lm(cr • 
gi), for all indices i. For the Grobner S-bases of P we have hence the following 
characterization. 

Proposition 3.8 (E-criterion). Let G be a Yj-basis of a T,-ideal I G P. Then, 
G is a Grobner Y,-basis of I if and only if for all f,g £ G,/,j / and for any 
cr, t £ E such that gcd(cr, t) — 1 and gcd(er • lm(/), r • lm(g)) =/= 1, the S-polynomial 
spoly(cr ■ f,T ■ g) has a Grobner representation with respect to E • G. 

Proof. Recall that G is a Grobner E-basis if and only if E • G is a Grobner basis 
of I. By the Buchberger's criterion [4] or by the Bergman's diamond lemma [2] 
this happens if and only if the S-polynomials spoly(er • /, r • g) have a Grobner 
representation with respect to E • G, for all /, g £ G, f, g ^ and cr, t £ E. By the 
product criterion we may restrict ourselves to considering only S-polynomials such 
that gcd(cr-lm(/),rdm(c/)) ^ 1 since lm(cr-/) — a -\xn(f) ,\vh(t • g) = T-lm(g). Then, 
let spoly(cr ■ f,T • g) be any such S-polynomial and put 5 = gcd(cr, r) and therefore 
o = 5a' , t = 5t' with a' , t' £ E, gcd(er', r') = 1. One has that spoly(er ■ /,r • g) = 
5 ■ spoly(cr' • /, t' ■ g) owing to Proposition EUl Note now that if spoly(er' ■ /, t' ■ g) = 
h = J2v fviy ■ gv) [y G E,/„ G P,g v G G) is a Grobner representation with 
respect to E • G then also spoly(tr ■ /, r • g) = 8 ■ h = J2 U (5 ■ fv){8v ■ g v ) is a Grobner 
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representation because -< is a E-ordering of P. We conclude that the S-polynomials 
to be checked for Grobner representations may be restricted to the ones satisfying 
both the conditions gcd(<r ■ lm(/), r • lni(g)) ^ 1 and gcd(cr, r) = 1. D 

The above result implies immediately the correctness of the following algorithm 
for the computation of Grobner E-bases. 

Algorithm 3.1 SigmaGBasis 

Input: H , a E-basis of a E-ideal I C P. 
Output: G, a Grobner E-basis of I. 
G:=H; 

B:={{f,g)\f,g€G}; 
while B / do 
choose (f,g) e B; 
B:=B\{(f,g)}; 

for all cr, t <G E such that gcd(er, r) = 1, gcd(<7 • lm(/), r • lva(g)) ^1 do 
h := REDUCE(spoly(cr • /, r • g), E • G); 
if /i ^ then 

B:=BU{(j,/ l ) 1 (/i,/»)|j6G}i 
G:=GU{/i}; 
end if; 
end for; 
end while; 
return G. 

Note that the function Reduce is given by the following standard routine. 



Algorithm 3.2 Reduce 



Input: G C P and / e P. 

Output: he P such that / -he (G) and h = or lm(/i) ^ LM(G). 

while /i ^ and lm(/i) e LM(G) do 

choose g e G,g ^ such that hn(g) divides lm(h); 

h:=h- (lt(h)/lt(g))g; 
end while; 
return h. 

Clearly, the chain criterion can be added to SigmaGBasis to shorten the number 
of S-polynomials that have to be reduced. In fact, one can view this algorithm as a 
variant of the Buchberger's procedure applied to the basis E • H, where Proposition 
3.81 provides the additional criterion "gcd(cr, r) = 1" to avoid useless pairs. Note 
also that for any pair of elements f,geG and for all a, r e E there are only a finite 
number of S-polynomials spoly(cr • /, r • g) satisfying both the criteria gcd(er, r) = 1 
and gcd(er • lm(/),r • lm(<?)) ^ 1. A proof is given by the arguments contained in 
Proposition ^. 71 

It is important to remark that we do not have general termination for the algo- 
rithm SigmaGBasis owing to Non-Noetherianity of the ring P. In fact, even if a 
E-ideal I C P has a finite E-basis this may be not true for its initial E-ideal LM(J) 



NOETHERIAN QUOTIENTS OF THE ALGEBRA OF PARTIAL DIFFERENCE ... 7 

that is I may have no finite Grobner S-basis. In the following sections we propose 
then two possible solutions to this problem. First, we introduce a suitable grading 
on P that is compatible with the action of £ which implies that the truncated 
version of the algorithm SigmaGBasis with homogeneous input stops in a finite 
number of steps. Another approach consists in obtaining finite Grobner E-bases 
when suitable elements belong to the given S-idcal /. More precisely, we prove the 
Noetherian property for suitable (quotient) E-algebras P/I. 

4. Grading and truncation 

A useful grading for the free S-algebra P can be introduced in the following way. 
Consider the set N = NU {—00} endowed with the binary operations max and +. 
Then (N, max, +) is clearly a commutative idempotent semiring. Moreover, for any 
°" = Hi a T G S we put deg(cr) = J2i a i- 

Definition 4.1. Let ord : M — >■ N be the unique mapping such that 
(i) ord(l) = —00; 

(ii) ord(r7jn) = max(ord(m), ord(n)), for any m,n£z M; 
(hi) ord(xi(o~)) = deg(cr), for all i > and a G S. 

Then, the map ord is a monoid homomorphism from (M, •) to (N, max). We call 
ord the order function of P. 

For any monomial m = a^^i)" 1 • • ■Xi k (Sk) ak different from 1 we have that 
ord(ra) = max(deg(#i), . . . ,deg(<5fc)). If Pd — ( m € M | ord(m) = d)x C P then 
P = © d6 ^ Pd is clearly a grading of the algebra P over the commutative monoid 

(N, max). 

Proposition 4.2. The following properties hold for the order function: 
(i) ord(er • m) = deg(cr) + ord(m), for any a 6 £ and m £ M ; 
(ii) ord(lcm(r7i, n)) — ord(mn) = max(ord(r7i), ord(n)), for allm,n G M. There- 
fore, if m I n then ord(m) < ord(n) . 

Proof. If m — 1 then ord(cr • m) = ord(m) = —00 = dcg(cr) + ord(m). If other- 
wise m — x il (Si) ai ■ ■ ■ Xi k (<5fe) Qfc then a ■ m = x il (oSi) ai ■ ■ ■ x ik (aSk) ak and hence 
ord(cr • m) = max(deg(<7(5i), . . . , deg(crS k )) = deg(cr) + max(deg(<5i), . . . , deg(4)) = 
deg(er) + ord(m). To prove (ii) it is sufficient to note that the order of a monomial 
does not depend on the exponents of the variables occurring in it. □ 

Definition 4.3. An ideal I c P is called ord-graded if I — ^2 d Id with Id — inPd- 
If I is also a Ti-ideal then a ■ Id C ^dcg(<r)+d f or an V c G S and rfgN. 

Let /, g G P, f ^ g he any pair of ord-homogeneous elements. Then, the S- 
polynomial h = spoly(/, g) is also ord-homogeneous and we have that ord(/i) = 
max(ord(/),ord(g)). If ord(/), ord(g) < d for some d G N, one has therefore that 
ord(/i) < d which implies the following result. 

Proposition 4.4 (Termination by truncation). Let I C P be a ord-graded H-ideal 
and let d £ N. Assume there is a ord-homogeneous basis H <Z I such that Hd = 
{/ £ H I ord(/) < d} is a finite set. Then, there exists also a ord-homogeneous 
Grobner Y,-basis G of I such that Gd is a finite set. In other words, if one uses for 
SigmaGBasis a selection strategy of the S-polynomials based on their orders then 
the d-truncated version of SigmaGBasis terminates in a finite number of steps. 
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Proof. In the algorithm SigmaGBasis one computes a subset G of a Grobner 
basis G' = £ • G obtained by applying the Buchberger's algorithm to the basis 
H' = £ • H of the ideal I. Moreover, Proposition 14.21 implies that the set H' and 
hence G' consists of ord-homogeneous elements. Define H' d = {a ■ f \ a G S, / G 
H, deg(cr) + ord(/) < d}. Since X^ = {cr G £ | deg(er) < d} and by hypothesis 
P/d are finite sets, we have that H' d C £d ■ Hd is also a finite set. Denote now 
Yd the finite set of variables of P occurring in the elements of H' d and define the 
subalgebra Pi d \ = K[Y d ] C P. In fact, the d-truncated algorithm SigmaGBasis 
computes a subset of a Grobner basis of the ideal Lj\ C Pm) generated by P/^. 
The Noetherianity of the finitely generated polynomial ring Pi d \ provides then 
termination. D 

Note that this result implies algorithmic solution to the ideal membership for 
finitely generated ord-graded E-ideals. Another consequence of the grading defined 
by the order function is that one has a criterion, also in the non-graded case, 
for verifying that a S-basis computed by the algorithm SigmaGBasis using a 
finite number of variables of P is a complete finite Grobner S-basis, whenever this 
basis exists. This is of course important because actual computations can be only 
performed over a finite number of variables. 

Definition 4.5. Let ~< be a monomial Y,-ordering of P. We say that -< is com- 
patible with the order function if ord(m) < ord(n) implies that m -< n, for all 
m,n G M. 

Proposition 4.6. Denote by -< the monomial H-ordering of P defined in Propo- 
sition [373\ and let < be the monomial ordering of Q — K[o~i, . . . , oy] that is used to 
define -<• If < is compatible with function deg then -< is compatible with function 
ord. 

Proof. Let m — m(5i) ■ • ■m{5k),n = n{5\) • • ■ n(6k) two monomials of P with 
m(Si),n(Si) G M(6i) and Si > . . . > S k (hence deg(<$i) > ... > deg(<Ji)). As- 
sume m -< n that is there is f < i < k such that m(Sj) = n(Sj) when j < i and 
m(5i) -< n(6i). If i > 1 or m(6i) ^ 1 one has clearly ord(m) = ord(n) = deg(5i). 
Otherwise, we conclude that ord(m) < deg(Si) = ord(n). D 

As before, we denote S^ = {a G S dcg(a) < d}. 

Proposition 4.7 (Finite S-criterion). Assume P be endowed with a monomial S- 
ordering compatible with the order function. Let G C P be a finite set and define the 
Ti-ideal I = (G)s- Moreover, denote d — max{ord(lm((?)) | g G G,g ^ 0}. Then, G 
is a Grobner T,-basis of I if and only if for all f,g G G and for any cr, r G £ such 
thatgcd(a,r) = 1 and gcd{a ■\m{f) 1 T ■\m.{g)) ^ I, the S-polynomial spoly(cr • f ,t ■ g) 
has a Grobner representation with respect to the finite set T,2d • G. 

Proof. Let spoly(cr • /, r • g) = h = ^2 V f v (y ■ g v ) be a Grobner representation with 
respect to E • G that is lm(/i) y lm(f u )(v ■ \m(g„)) for all v. We want to bound the 
degree of the elements i/£S occurring in this representation. Put m = lm(/), n = 
lm(<?) and hence lm(cr • /) = a -m, lm(cr • g) = a -n. By product criterion, we assume 
that u = gcd(er-m, r-n) ^ 1. Then, there is a variable Xi(aa) = Xi(r/3) that divides 
u where Xi{a) divides m and hence deg(a) < ord(m) < d and Xi{j3) divides n and 
therefore deg(/3) < ord(n) < d. Then act = t/3 and one has that a \ (3, r \ a because 
we assume also that gcd(tr, r) = 1. We conclude that deg(cr), deg(r) < d and if v = 
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lcm(<7 • m, t ■ m) then ord(i>) = max(dcg(er)+ord(77j),deg(T)-|-ord(n)) < 2d. Clearly 
v >- lm(/i) y i/-lm(g v ) and hence 2d > ord(w) > deg(t , )+ord(lm(g I/ )) > dcg(^). D 

Under the assumption of a S-ordcring compatible with the order function and 
for S-idcals that admit finite Grobner E-bases, by above criterion one has an actual 
algorithm to compute such a basis in a finite number of steps. In fact, this can be 
obtained as an adaptative procedure that keeps the bound 2d for the degree of the 
elements of £ applied to the generators, constantly updated with respect to the 
maximal order d of the leading monomials of the current generators. 

5. AN ILLUSTRATIVE EXAMPLE 

In this section we apply the algorithm SigmaGBasis to an example arising from 
the discretization of a well-known system of partial differential equations. Consider 
the unsteady two-dimensional motion of an incompressible viscous liquid of constant 
viscosity which is governed by the following system 

U X + Vy = 0, 

U t + UU X + VUy +P X - j{U XX + Uyy) = 0, 

Vt + UV X + Wy +Py - ^(V XX + Vyy) = 0. 

The last two nonlinear equations are the Navier-Stokes equations and the first linear 
equation is the continuity one. Equations are given in the dimensionless form where 
(u, v) represents the velocity field and the function p is the pressure. The parameter 
p denotes the Reynolds number. For defining a finite difference approximation of 
this system one has therefore to fix X = {u,v,p} and £ = (<7i,<72, U3) since all 
functions are trivariate ones. To simplify notations, we identify £ with the addi- 
tive monoid N 3 and we have that P = K[X(E)] = K[u(i,j,k),v(i,j,k),p(i,j,k) | 
hj, k > 0]. The base field K is the field of rational numbers. The approximation 
of the derivatives of the function u is given by the following formulas (forward 
differences) 

^ u{x+h x .y.t)-u{x.y.i) _ m(1.0.0)-m(0.0.0) 
Ux ~ h* ~~ h x 



II, 



u( X ,y+h y ,t)-u( X ,y,t) _ m(0,1,0)-m(0,0,0) 

V ~ hy - hy ' 

u(x,y,t+h t )-u(x,y,t) _ m(0,0,1)-m(0,0,0) 



U t h t h t ' 

^ u(x+2hy,y,t)-2u{x+h x .y.t)+u(x,y,t) _ m(2,0,0)-2m(1,0,0)+m(0,0,0) 
U XX J^ — J^S , 

^ u(x,y+2h v ,t)-2u(x,y+hy,t)+u(x,y,t) _ m(0.2,0)-2m(0,1.0)+m(0.0,0) 

Uyy p^ — jp , 

y y 

where h x ,h y ,ht are parameters (mesh steps). One has similar approximations for 
the derivatives of the functions v,p. Then, the Navier-Stokes system is approxi- 
mated by the following system of partial difference equations 

/i = f (u(l,0,0)-u(0,0,0)) + f (v(0,l,0)-v(0,0,0)) = 0, 

h = -^r«(2, 0, 0) - ^«(0, 2, OH (Jj + £«(0, 0, 0))«(1, 0, 0) + 

£p(l, 0, 0) + (^ + £ v(0, 0, 0))u(0, 1, 0) + £u(0, 0, 1) - £u(0, 0, 0) 2 - 

(^ + ^ + 717 + £v(P> ^ 0))«(0, 0, 0) - £p(0, 0, 0) = 0, 

h = ~^v(2, 0, 0) - ^v(Q, 2, 0) + (^ + £«(0, 0, 0))v(l, 0, 0) + 

(w + t^ (0 ' °' 0))w(0 ' *> 0) + t p(0 ' X ' °) + ^ w(0 ' °' X) " ^ m(0 ' °' 0)w(0 ' °' 0) 
[ ~£v(0, 0, 0) 2 - (^ + ^ + £*)„((), 0, 0) - ^(0, 0, 0) = 0. 
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We encode this system as the E-ideal I = {fi, fi, h)n C P and we want to com- 
pute a (hopefully finite) Grobner E-basis of /. For instance, we may want to have 
such basis to check for "strong-consistency" [10J of the finite difference approxima- 
tion. We fix the degree reverse lexicographic ordering on the polynomial algebra 
K[a\,cr<2,a^\ (ai > a 2 > 173) and the lexicographic ordering on K[u, v,p] (u >- v >- 
p). By Proposition 13.31 one obtains therefore a (block) monomial E-ordering for P. 
Note that this ordering is compatible with the order function and hence Proposition 
14.71 is applicable to certify completeness of a Grobner E-basis computed over some 
finite set of variables {u(i,j,k),v(i,j,k),p(i,j,k) | i+j + k < d}. 

With respect to the monomial ordering assigned to P, the leading monomials of 
the ^-generators of / are lm(/i) = u(l, 0, 0),lm(/ 2 ) = u(2, 0, 0),lm(/ 3 ) = v(2, 0, 0). 
Since o\ ■ lm(/i) = lm(/2), by interreducing fi with respect to the set E • {/1, fa} 
we obtain the element 

ft = h v htv(l, 1,0) - h x h t u(0, 2, 0) - hyhtv(l, 0, 0) + ph 2 y h tP (l, 0, 0)+ 
(2h x h t + ph x h y h t v(0, 0, 0))ti(0, 1,0) - (h v h t + ph x hyh t u(0, 0, 0))«(0, 1,0)+ 
P h x h 2 y u{Q, 0, 1) - (h x h t + ph x hl)u(0, 0, 0) - ph 2 y h t p{Q, 0, 0) + h y h t v{0, 0, 0) 

whose leading monomial is lm^) = v(l, 1, 0). Owing to the E-criterion, the only 
S-polynomial to consider is then spoly(eri • f 2 , o 2 ■ Is) whose reduction with respect 
to E • {/1, f 2 , fs} leads to the new element 

U = h 3 y h t p{2, 0, 0) + (-h 2 x h y htu(0, 0, 0) - h 3 x h t v(0, 0, 0) + h x htv(0, 1,0)+ 
h 2 x h y h t u{0, 1, 0))u(0, 2, 0) + {-hlh y h t v(l, 0, 0) + h x h v h t v(0, 1, 0))u(0, 2, 0) + 
h 2 x h y h tP (0, 2, 0) + (2h 2 x h y h t v(0, 1, 0) + 2h x h 2 y h t u{0, 1, 0) - h x h v h t v(0, 0, 0)- 
2h x h 2 h t u(0, 0, 0))u(l, 0, 0) + (phlhlhtv(0, 0, 0) + ph x fVyh t u(0, 0, 0)- 
ph 2 x h 2 h t v{0, 1, 0) - ph x h 3 y h t u(0, 1, 0) - 2h 3 y h t )p(l, 0, 0) - (ph x h? y h t v(0, 0, 0)+ 
2h 2 x h y h t )u(0, 1, 0) 2 + ((ph 2 x h 2 y h t u(0, 0, 0) - ph 3 x h y h t v(0, 0, 0) - 2h 3 x h t )v(0, 1, 0)- 
ph 2 x hlu(0, 0, 1) + (3h 2 x h y h t + ph 2 x h 3 y + ph 2 x h 2 y h t v{Q, 0, 0))u(0, 0, 0)+ 
ph 3 x h y h t v(0, 0, 0) 2 + (2/i^ _ 2h x h 2 y h t )v(0, 0, 0) + ph x h 3 y h tP (0i 0, 0))«(0, 1, 0)+ 
phlh y htu{0, 0, 0)u(0, 1, 0) 2 + (-p/i3/i 2 u (0, 0, 1) - ph x h 2 y htu{Q, 0, 0) 2 + 
(ph 3 x h 2 y + /i^/it - ph 3 x h y h t v(0, 0, 0))u(0, 0, 0) - 3h 2 x h y h t v(0, 0, 0) + 
ph 2 x h 2 y h tP (0 : 0, 0))«(0, 1, 0) - 2h 2 x h y htp(0 : 1, 0) + (pft*^ w ( 0) 0, 0)+ 
phlhlulo, 0, 0))u(0, 0, 1) - (p/i 2 ^ + h 2 x h y h t )u(0, 0, 0) 2 + ((-/i^/i* - ph 3 x h 2 y + 
2h x h 2 h t )v(0, 0, 0) - ph x h 3 y h t p(0, 0, 0))u(0, 0, 0) + 2h t h 2 x h y h t v(0, 0, 0) 2 - 
ph 2 x h 2 h t v(0, 0, 0)p(0, 0, 0) + (h 3 y ht + h 2 x hyht)p(0, 0, 0). 

The leading monomial of this difference polynomial is lm(/4) = p(2,0, 0) and no 
more S-polynomials are defined. We conclude that the set {/1, f' 2 , ^3, f^} is a (finite) 
Grobner E-basis of the E-ideal I C P. Since we make use of a monomial E-ordering 
for P, this is equivalent to say that E • {/1, f 2 ,h, ft} is a Grobner basis of the ideal 
/ and this can be verified also by applying the usual Grobner bases routines to the 
basis E • {fi, fz, fz\- Because the maximal order in the input generators is 2, by 
Proposition 14.71 it is reasonable to bound initially the order of the variables of P 
to 4 or 5. Note that the computing time for obtaining a Grobner basis of / with 
implementation in Maple of the Faugere's F4 algorithm amounts to 20 seconds for 
order 4 and 5 hours for order 5 on our laptop Intel Core 2 Duo at 2.10 GHz with 
8 GB RAM. By the algorithm SigmaGBasis that we implemented in the Maple 
language as a variant of the Buchberger's algorithm (see [T5]), the computing time 
for a Grobner E-basis of / is instead seconds for order 4 and 3 seconds for order 
5 since just two reductions are needed. 
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6. A NOETHERIANITY CRITERION 

As already noted, a critical feature of the algebra of partial difference poly- 
nomials P = K[X(E)] is that some of its E-ideals are not only infinitely gen- 
erated but also infinitely ^-generated. One finds an immediate counterexam- 
ple for E = (a) that we identify with the additive monoid N. Then, the ideal 
/ = (xi(0)Xi(l),Xi(Q)xi(2), . ..)s has clearly no finite N-basis. It is interesting 
to note that if one considers a larger semigroup like Inc(N) = {/ : N — >• N | 
/ strictly increasing} acting on P = K[X{H)] as / • Xi(j) = Xi(f(j)), one obtains 
that P is Inc(N)-Noetherian [T. We may say hence that the monoid N is "too small" 
to provide N-Noetherianity. One way to solve this problem is to consider suitable 
quotients of the algebra of partial difference polynomials where Noetherianity is 
restored. A similar approach is used for the free associative algebra which is also 
Non-Noetherian where the concepts of "algebras of solvable type, PBW algebras, 
G-algebras", etc naturally arise (see for instance [2fJ ). 

We start now with a general discussion for (commutative) algebras generated by 
a countable set of elements. Let Y — {2/1,2/2, • ■ •} be a countable set and denote 
P = K[Y] the polynomial algebra with variables set Y. Since P is a free algebra, all 
algebras generated by a countable set of elements are clearly isomorphic to quotients 
P' = P/J, where J is some ideal of P. To control the cosets in P', a standard 
approach consists in defining a normal form modulo J associated to a monomial 
ordering of P. Then, let -< be a monomial ordering of P such that j/i -< 2/2 -^ 

Definition 6.1. Put M = Mon(P) and denote M" = M\\m(J). Moreover, define 
the subspace P" = (M")k C P. The elements of M" are called normal monomials 
modulo J (with respect to -<). The polynomials in P" are said in normal form 
modulo J. 

Since P is endowed with a monomial ordering, by a standard argument based 
on the procedure Reduce applied for the set J one obtains the following result. 

Proposition 6.2. A K-linear basis of the algebra P' is given by the set M' — 
{m + J\me M"}. 

Definition 6.3. Let f E P. Denote NF(/) the unique element of P" such that 
f - NF(/) e J. In other words, one has NF(/) = ReduCe(/, J). We call NF(/) 
the normal form of / modulo J (with respect to -<). 

By Proposition 16.21 one has that the mapping / + J 1— »• NF(f) defines a linear 
isomorphism between P' = P/J and P" — (M")k- An algebra structure is defined 
hence for P" by imposing that such mapping is also an algebra isomorphism that 
is we define f ■ g — NF(fg), for all /, g 6 P". By means of this, one has a complete 
identification of M', M" and P', P" that we will use from now on. We define then 
the set of normal variables 

Y' = Y C\M' = F\lm(J). 

Clearly, normal variables depend strictly on the monomial ordering one uses in P. 

Proposition 6.4 (Noetherianity criterion). Let P be endowed with a monomial 
ordering. If the set of normal variables Y' is finite then P' is a Noetherian algebra. 

Proof. It is sufficient to note that all normal monomials are product of normal 
variables and therefore the quotient algebra P' = P/J is in fact generated by the 
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set Y'. If Y' is finite then P' is a finitely generated (commutative) algebra and 
hence it satisfies the Noetherian property. □ 

We need now to introduce the notion of Grobner basis for the ideals of P' = P/J. 
After identification of cosets with normal forms, recall that M ' — M \ lm( J) and 
P' = (M')k is a subspace of P endowed with multiplication / • g = NF(fg), for 
all f,g £ P'. Then, all ideals V C P' have the form V = I/J = {NF(/) | / £ I}, 
for some ideal J £ I £ P. Note that NF(/) G / for any / £ I, which implies 
that in fact I' = I n P' . Since the quotient algebra P'/I' is isomorphic to P/I and 
Grobner bases give rise to X-linear bases of normal monomials for the quotients, 
one introduces the following definition. 

Definition 6.5. Let I' = I n P' be an ideal of P' where I D J is an ideal of P. 
Moreover, let G' C V . We call G' a Grobner basis of I' if G' U J is a Grobner 
basis of I . 

Let G c P. Recall that LM(G) denotes the ideal of P generated by the set 
lm(G) = {lm(g)\g£G,g^O}. 

Proposition 6.6. Let V be an ideal of P' and let G' C I'. Then, the set G' is a 
Grdbner basis of I' if and only z/LM(G") = LM(I'). 

Proof. Assume LM(G') = LM(I'). Let f £ I and denote /' = NF(/). If lm(/) £ 
LM(J) then lm(/) = lm(/'). Since lm(/') £ LM(J') C LM(G') then lm(/) = 
lm(/') = mlm(g'), for some m £ M,g' £ G'. We conclude that G' U J is a 
Grobner basis of /. Suppose now that the latter condition holds. Since G' C I' 
we have clearly that LM(G') C LM(I'). Let now /' £ V C /. Then, there is 
m £ M,g' £ G' U J such that lm(/') = mlm(g). Since lm(/') £ M' then also 
lm(fir) £ M' and hence g £ G '. □ 

Proposition 6.7. Assume that the set of normal variables Y' = Y n Af ' is finite. 
Then, any monomial ideal I = (I P\ M') C P /ias a finite basis. 

Proof. It is sufficient to invoke the Dickson's Lemma (see for instance [5]) for the 
ideal I which is generated by normal monomials that are products of a finite number 
of normal variables. □ 

Corollary 6.8. If Y' is a finite set then any ideal I' C P' has a finite Grdbner 
basis. 

Proof. According to Proposition ^. 61 consider the ideal LM(I') C P which is gener- 
ated by the set of normal monomials lm(P). Then, it is sufficient to apply Propo- 
sition [23 to this ideal. □ 

It is clear that if G is any Grobner basis of an ideal J ^ P then Y' = Y\ lm(G). 
Since Y is a countable set, if Y' is a finite set and hence P' = K\Y'] is a Noetherian 
algebra then G needs to be an infinite set. In general, such Grobner basis cannot 
be computed but this is possible when P is a S-algebra and J is a S-ideal. From 
now on, we assume again that P = K [X(S)] is the algebra of partial difference 
polynomials. Let J C P be any E-ideal that is P' = P/J is any S-algebra which 
is (finitely) ^-generated by -X"(l) and J is the ideal of E-algebra relations satisfied 
by such generating set. Let P be endowed with a monomial E-ordering -< and 
define, as before, the set M' C M = Mon(P) of all normal monomials and the set 
A(S)' = A(S) n M' of all normal variables. Then P' is an algebra generated by 
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X(T,) f . More precisely, the S-algebra P' is E-generated by X(l)' = X(l) n M' . In 
fact, the set X(E) \ X(E)' = -X"(S) fl lm(J) is clearly stable under the action of S 
that is if Xi(a) £ X(T,)' then ^(1) £ -X"(l)'. The structure of E-algebra of P' is 
given by defining, for all a £ E and 2^(1) € -X"(l)' 

ff-ai(l)=NF(u(ff)). 

Proposition 6.9 (Finiteness criterion). TTie sei of normal variables X(T,)' is finite 
if and only if for all 1 < i < n, 1 < j < r one /ias t/iat Xi(cr lJ ) £ lm(J), /or some 
integers d^ > 0. 

Proo/. Put ^(E) = {^(cr) I a £ E} and denote ^(E)' = ^(E) n X(E)', for any 
i = 1, 2, . . . , n. We have then to characterize when Xj(E)' is a finite set. Consider 
the polynomial algebra Q — K[a\ 1 . . . , ay] and a monomial ideal I C Q. It is well- 
known that the quotient algebra Q / 1 is finite dimensional if and only if there are 

integers dj > such that a J £ /, for all j = 1, 2, . . . ,r. In other words, £,(E)' is a 

finite set if and only if there exist integers dij > such that Xi(a^' 3 ) £ lm(J), for 
all indices j, j. □ 

Corollary 6.10 (Termination by membership). Let J c P be a Yi-ideal such that 
Xiia, 1 ' 3 ) £ lm(J) for all 1 < i < n. 1 < j < r and for some dij > 0. Then J has a 
finite Grobner "E-basis. 

Proof. Denote I = {x l {af 1 ) \ 1 < i < n, 1 < j < r) E and L = LM(J). Then, we 
have that I C L and the ideal P/I C P/I has a finite basis owing to Proposition 
16 .71 and Proposition 16. 91 In other words, the E-ideal L has a finite E-basis given by 
the finite E-basis of I together with the finite basis of L/I. □ 

Note that the above result is generally not a necessary condition for finiteness of 
Grobner S-bases. Consider for instance the example presented in Section 5 of |18j . 
Nevertheless, Corollarv l6.10l guarantees termination of the algorithm SigmaGBasis 
when a complete set of variables Xj(er-* 3 ) for all i,j, occur as leading monomials of 
some elements of the Grobner E-basis at some intermediate step of the computation. 
Of course, if such elements belongs to the input E-basis of a E-ideal J C P one 
knows in advance that all properties of Noetherianity and termination are provided 
for the quotient P' = P/J. In particular, one may have that the polynomials 
fij £ P such that ]m(fij) = Xi(o- IJ ) are themselves a Grobner E-basis of J. This 

is clearly the case when J = (xi(<7* J ) \ 1 < i < n,l < j < r)s- For all d > 0, define 
now 

J (d) = (xi(a) | 1 < i < n,deg(er) = d + l) s D {x l {aj +l ) \l<i<n,l<3< ^}s 

and put j(-°°) = (X(l)) s = (X(E)). If P = ® deS ,P d is the grading of P defined 
by the order function then the subalgebra P^ = ® i<d Pi C P is clearly isomorphic 
to the quotient P j 3d and hence it can be endowed with the structure of E-algebra. 
In other words, to make use of the Noetherian (E-)subalgebras filtration 

K = p(-oo) ,_ p(0) c p(l) c c p 

to perform concrete computations with Grobner E-bases as explained in Section 4 
corresponds to work progressively modulo the E-ideals 

(X(T)) - J*-- 00 '* D J {0) D J (1) D...D0 
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providing a finite number of normal variables. 

Another relevant case is the ordinary one that is when £ = (a). In this case, 
any set of polynomials fi, ■ ■ ■ , f n € P such that lm(/j) = Xi{(j di ) (di > 0) is a 
Grobner S-basis since no S-polynomials are defined (up to product criterion). One 
has finally the following interesting case. 

Proposition 6.11. Assume that K is a field of constants and consider the linear 
polynomials fy = J2o<k<d i:j CijkXi{o~j) e P where c ijk e K and c l3dl] = 1, for all 

1 < i < n, 1 < j < r. Then lm(/,y) = Xi(o~ - %] ) and the set {fij} is a Grobner 
Yi-basis. 

Proof. Since X(E) is endowed with a S-ordering, one has that Xi{a\j) -< Xi(a l A if 

k < I and hence lm(/jj) = Xi{a^ J ). Then, the only S-polynomials to be considered 
are 

a = spo\y{af"-f ip ,af'-Si q )= ^ ^pkMrf" a p) ~ X! c iq iXi{af"<j q ), 

Q<k<d ip 0<Kd iq 

r k 



for all 1 < i < n and I < p < q < r. By reducing s with polynomials cc • fi q and 

q ' tip 



o~ ■ fip one obtains 



/ J Ci p kCi q iXi(a q a p ) + 22 Ci q ic ip kXi(o-pO- q ) = 0. 

0<k<d ip fi<l<d iq Q<l<d iq ,0<k<d ip 



Note explicitely that if £ acts on K in a non-trivial way then generally 
s' = - Yl ^t" ' c ipk)(°l ■ ^qi)^{v l q ° k p ) + 

0<k<d ip ,0<Kd 



D 



p) u^-i\u,iq 



Yl ( a P P ■ c ni)^\ ■ Cipk)Xi{a k p(j q ) + 0. 



0<l<d ia ,0<k<di. 



7. A NOETHERIAN S-ALGEBRA OF SPECIAL INTEREST 

From now on we assume that K is a field of constants. We define the ideal 
J = (fij)-z C P where fy = Y,o<k<d tj c ijkXi{a^) {c ijk G K,c ljdi] = 1), for any 
1 < i < n,l < j < r. We want to describe the (Noetherian) E-algebra P' = P/J. 
To simplify notations and since they are interesting in itself, we consider separately 
the cases when r = 1 and n = 1. Then, assume S = (cr) and hence P' = P/J 
where J = (/i, . . . ,/„)s with /; = I] <fe<<i, c Jfe^( cr ' c ) ( c ife € K,dd t = 1). Define 
Q = K[o~] the algebra of polynomials in the single variable a and denote gi — 
Tlo<k<d- c ik ak ^ Q- Moreover, put d — '^ li di and let V = K d . Finally, consider 
the d x d block-diagonal matrix 



A = A 1 ®...(BA n = 



( A 1 ... \ 

A 2 ... 



V ... A n J 
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where each block Ai is the companion matrix of the polynomial gi that is 

/ ... -cm \ 
1 ... -cii 
1 



Ai = 



-Cii 



V ... 1 -Ctf-1 / 

Note that A has all entries in K and it can be considered as the Frobenius normal 
form of any d x d matrix provided that gi | ... | g n . Then, the monoid £ or 
equivalently the algebra Q acts linearly over the vector space V by means of the 



representation a 



(a k ) 



where q 



A . If {w 9 }i< 9 <d is the canonical basis of V, we denote 
= Ylj<i dj + k + 1 for all 1 < i < n, < k < d t . We have 
hence Xi(a k ) = A k Xi(l) = a k ■ xi(l). In other words, for the Q-module V one 
has the decomposition V = ® ■ Vi where Vi is the cyclic submodule generated by 
Xi(l) and annihilated by the ideal (gi) C Q. Denote now by R' the (Noetherian) 
polynomial algebra generated by the finite set of variables X(S)' = {xi(a k ) | 1 < 
i < n,0 < k < d{\ that is V coincides with the subspace of linear forms of R' . 
Then, one extends the action of the monoid S = (a) to the polynomial algebra R' 
in the natural way that is by putting, for all k > and Xi(a°) 6 X(S)' 

a k -x i (a j )=A k x i (a j ). 

Denote by end(P) the algebra of all X-linear maps P — > P. Note that the 
representation p : S — > End(P) can be extended linearly to p : Q — > end(P). Then, 
one has that fc = Y.k c ik x i{ ak ) = J2k °i^ k -Xi(l) = ^-^(1), for all i = 1, 2, . . . ,n. 

Proposition 7.1. IfS = (a) then the Y,-algebras P',R' are ^.-isomorphic. 



Proof. By Proposition 16. 1 II we have that the set {/»} is a Grobner E-basis of the 
E-ideal J C P and it is clear that the set of normal variables modulo J is exactly 
X(T,)' = {x t (o- k ) \ 1 < i < n,0 < k < d l }. Moreover, since R' C P and /* = 
g % ■ Xi(l) one has that NF(a; 4 (cr fe )) = NF(cr fc • x z (l)) = A^'x^l), for all k > and 
^(1) G X(l)'. D 

Note that R' is S-generated by the set X(l)' = {xi(l) \ 1 < i < n,di > 0}. 
Since P is a free E-algebra, a surjective E-algebra homomorphism tp : P — >• i?' is 
defined such that 

x»(l) if dj > 0, 
otherwise. 



Xi(l) 



Then, the above result states that the E-ideal Ker ip C P oi all S-algebra relations 
satisfied by the generating set X(l)' U {0} of R' is exactly J. 

Let now S = (a\ , . . . , ay) and suppose X = {x}. Then P' = P/J where 
J = (fi,---,fr)z with fj = Y,a<k<d i c 3k x ( a 'j) ( c Jk G K , c jdj = 1). Define 
Q = K [(7i, . . . ,<7 r ] the algebra of polynomials in the variables <7j and denote 
<?j = X)o<fc<d C jfc (7 j : £ Q- One has clearly that /j = g^ • x(l). As before, we 
consider the companion matrix Aj of the polynomial gj in the single variable o~j . 



If d = Y[ j dj then the monoid S = Ei 



Q = Ql®---®Qr(Qj = K[<Tj 

the representation 



x E r (Ej = (aj)) that is the algebra 
K d my means of 



acts linearly over the space V 

a kr ^ A kl ®---®A kr 7 
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where A^ 1 • • • <g> A kr denotes the Kronecker product of the matrices A- 3 . In other 
words, the Q-module V is the tensor product V — V\ <E) ■ ■ ■ <8> V r where Vj is the 
cyclic Qj-module defined by the representation a k (-> A k j. If {wfe 1 ® • • ■®Vk r }i<k i <d :j 
is the canonical basis of V, we put x(a 1 1 ■ ■ -o- kr ) = Vk 1 +i ® ■ ■ • ® Vfe r +i, for all 
1 < J < r , < fcj < dj . One has then 

x(cr kl ■ ■ • tr*') = (A* 1 ® . . . ® A^)a;(l) = (a^ 1 ■ • -a*-) • x(l), 

that is V is a cyclic module generated by a;(l). Denote now by R' the polynomial 
algebra generated by the finite set of variables -X"(E)' = {x(a l 1 ■ ■ ■ a kr ) | 1 < j < 
r, < kj < dj} that is V is the subspace of linear forms of R'. Again, we extend 
the action of the monoid E = (ci, • • • , ^r) to the polynomial algebra i?' by putting, 
for all ki,...,kr>0 and aj(<j) £ X(Y,)' 

(a kl ■ ■ ■ o*T) • x(a) = (^ ® • • • ® A^)x(<r). 

Proposition 7.2. If X = {x} then P',R' are E- isomorphic. 

Proof. Assume d ^ that is dj y^ for all j. Again, by Proposition 16.111 one has 
that the set {fj} is a Grbbner E-basis of J C P and the set of normal variables 
modulo J is clearly X(E)' = {^(crf 1 • • ■ a k /) \ 1 < j < r,0 < kj < dj}. Moreover, 
because R' C P and fj=9j-x(l) we obtain that, for all hi, . . ■ ,k r > 

NF(s(a^ • • • o*")) = NF((ct^ • • • a^) • a:(l)) = (A kl ® . . . ® A^)x(l). 

Finally, if d = then P' = R' = K. D 

Note that for d ^ one has that R' is E-generated by the element x(l). Then, 
the above result implies that the E-ideal J c P coincides with the ideal of E- 
algebra relations satisfied by the generator x(l) that is the kernel of the E-algebra 
epimorphism P — > R' such that x(l) H> x(l). 

Consider finally the general case for the E-algebra P' — P/J where J = (fij)s 
and fij = J2a<k<d tj c vkXi(o- k ) with c tjk e K,Ci jdii = 1, for all 1 < i < n and 
1 < 3 < r - By combining the previous results, it is clear that such structure arises 
from the Q-module V — K d where d = Yli<i<n Yii<i<r ^v an< ^ the representation 
is given by the mapping 

3 i 3 

where Aij is the companion matrix of the polynomial gij = J2o<k<d c ijkO~ k - ln 
other words, we have that V = ®^ (££) . Vij where Vij is the cyclic Qj-module anni- 
hilated by the ideal (g^) c Qj. By denoting Xi(l) the generator of the Q-module 
(g) Vij, we obtain that P' is isomorphic to the E-algebra R! = K[X(Y>)'] where 

X(E)' = {x l (a k ' 1 ■ ■ -o- kir ) \ 1 < i < n,l < j < r,0 < kij < d io } is the canonical 
basis of the space V. Then, one has that J = (fij)s is exactly the E-ideal of 
E-algebra relations satisfied by generating set X{1)' U {0} of R' . 

8. Another example 

With the aim of showing some application of the previous results, we fix now 
a setting that has been already considered in [25]. Note that in our approach all 
computations can be performed over any field. Fix r = 1 that is E = (a) and 
Q = K[a]. Consider S^ the symmetric group on d elements and let 7 E §d be any 



NOETHERIAN QUOTIENTS OF THE ALGEBRA OF PARTIAL DIFFERENCE . . . 
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permutation. Denote T = (7) C Sd the cyclic subgroup generated by 7. Moreover, 
let 7 = 71 • • • 7„ be the cycle decomposition of 7 and denote di the length of the cycle 
7j. Consider the polynomial algebra R 1 — i^x^c- 7 ) | 1 < i < n, < j < di] and 
identify the subset {xi(l), . . . , Xi(a di ~ 1 )} with the support of the cycle %. Clearly 
R' is a T-algebra that is there is a (faithful) group representation p' : T — > Aut(R') 
Consider now the polynomials gi 
matrix 



a di — 1 e Q and define the dx d block-diagonal 



A = 



Mi 






A 2 



V ... A n J 

where each block Ai is the companion matrix of the polynomial gi that is the 
permutation matrix 

/ ... 1 \ 
1 ... 
A . = 1 ... 



,£„(!),..., X„(cr d " X ) 



\ ... 1 / 

If we order the variables of R' as xi(l), . . . ,xi(a dl ^ 1 ), . 
then the representation p' is defined as -f k ■ Xi(a^) — A k Xi(a : '), for all i,j,k. 
In other words, one has that R 1 is a E-algebra isomorphic to P' = P/J where 
J = (A,...,/„)s and fi = g t ■ Xi(l) = x t (a dt ) - Xi(l) e P. Consider now a 
T-ideal (equivalently a E-ideal) L' = (hi,...,h m }r C R 1 and define the E-ideal 
L = (hi, . . . , h m , /1, . . . , / n )s C P. According with Definition 16.51 and the iden- 
tification of R' with the quotient P' one has that G' C £' is a Grobner T-basis 
(equivalently E-basis) of £' if by definition G" U {/1, ...,/„} is a Grobner S-basis of 
i. In practice, the computation of G' is obtained by the algorithm SigmaGBasis 
that generalizes the Buchberger's procedure to the free S-algebra P. Note that 
such computation terminates owing to Corollary 16. 101 

To illustrate the method we fix now 7 = (12345678) £ §s an d K — Q. To simplify 
the variables notation we identify E with N that is R' = K[x(0),x(l), . . . ,x(7)]. 
Consider the following T-ideal of R' 

V = (x(0)x(2) - x(l) 2 ,x(0)x(3) - x(l)a;(2))r = 

(s(0)a;(2) - x(l) 2 , x(l)x(3) - x(lf , a;(2)x(4) - x(3) 2 , a;(3)x(5) - x(4) 2 , 

x(4)x(6) - a;(5) 2 , x(5)x(7) - x(6) 2 ,x(7) 2 - x(0)x(6),x(l)x(7) - x(0) 2 , 

a;(0)x(3) - x(l)a;(2),a;(l)x(4) - a;(2)o;(3),a;(2)x(5) - a;(3)x(4), 

x(3)x(6) - a;(4)o;(5),a;(4)x(7) - x(5)a;(6),a;(6)a;(7) - a;(0)x(5), 

a;(0)x(7) - x{\)x(^) 1 x(2)x{7) - aj(0)a;(l)>. 

Note that a;(0)a:(2) - a;(l) 2 ,a;(l)x(3) - a;(2) 2 ,a;(0)x(3) - x(X)x(2) are well-know 
equations of the twisted cubic in P 3 . Define now / = x(8) — x(0) 6 P and hence 
R! = P' = P/J where J = (/) E . Then, a Grobner T-basis (or E-basis) of V is 
obtained by computing a Grobner E-basis of the ideal 

L = (x(0)x(2) - a;(l) 2 , x(0)x(3) - z(l)x(2), /) E C P. 

Fix for instance the lexicographic monomial ordering on P (hence on R') with 
x(0) -< x(l) -< . . . which is clearly a E-ordering. The usual minimal Grobner basis 
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of L' consists of 54 elements whose leading monomials are 

x(7) 2 ,x(6)x(7), 

x(0)x(2) -> x(l)x(3) -> x(2)x(4) -> x(3)x(5) -> x(4)x(6) -> x(5)x(7), 

x(0)x(3) -> x(l)x(4) -► x(2)x(5) ->■ x(3)x(6) -> x(4)x(7),x(2)x(7), 

x(l)x(7),x(0)x(7),x(6) 3 ,x(0)x(4) 2 -> x(l)x(5) 2 -> x(2)x(6) 2 , 

x(0) 2 x(4) -> x(l) 2 x(5) -» x(2) 2 x(6) -> x(3) 2 x(7),x(0) 2 x(6), x(0)x(6) 2 , 

x(l)x(6) 2 ,x(l) 2 x(6),x(3) 2 x(4) -> x(4) 2 x(5) -> x(5) 2 x(6), 

x(4)x(5) 2 -^x(5)x(6) 2 ,x(0)x(1)x(6),x(0)x(4)x(5) -> x(l)x(5)x(6), 

x(0)x(5)x(6),x(l)x(2)x(6),x(2) 4 -> x(3) 4 -S- x(4) 4 -> x(5) 4 ,x(0) 3 x(5), 

x(0)x(5) 3 , x(2) 3 x(3), x(2)x(3) 3 ->• x(3)x(4) 3 , x(0) 2 x(5) 2 , x(0) 2 x(l)x(5), 

x(2) 2 x(3) 2 ,x(l) 2 x(2) 3 ,x(l) 4 x(2) 2 ,x(l) 6 x(2),x(l) 8 . 

Note that the arrow between two monomials means that a monomial can be ob- 
tained by the previous one by means of the E-action. Then, the minimal Grobner 
T-basis of L' has just 32 elements and their leading monomials are 

x(7) 2 , x(6)x(7), x(0)x(2), x(0)x(3), x(2)x(7), x(l)x(7), x(0)x(7), x(6) 3 , 

x(0)x(4) 2 , x(0) 2 x(4), x(0) 2 x(6), x(0)x(6) 2 , x(l)x(6) 2 , x(l) 2 x(6), x(3) 2 x(4), 

x(4)x(5) 2 , x(0)x(l)x(6), x(0)x(4)x(5), x(0)x(5)x(6), x(l)x(2)x(6), x(2) 4 , 

x(0) 3 x(5), x(0)x(5) 3 , x(2) 3 x(3)x(2)x(3) 3 , x(0) 2 x(5) 2 , x(0) 2 x(l)x(5), 

x(2) 2 x(3) 2 ,x(l) 2 x(2) 3 ,x(l) 4 x(2) 2 ,x(l) 6 x(2),x(l) 8 . 

Precisely, such basis is given by following elements 

x(7) 2 - x(0)x(6), x(6)x(7) - x(0)x(5), x(0)x(2) - x(l) 2 , x(0)x(3) - x(l)x(2), 

x(2)x(7) - x(0)x(l), x(l)x(7) - x(0) 2 , x(0)x(7) - x(l)x(6), x(6) 3 - x(0)x(5) 2 , 

x(0)x(4) 2 - x(2)x(3) 2 , x(0) 2 x(4) - x(l) 2 x(2), x(0) 2 x(6) - x(0)x(l)x(5), 

x(0)x(6) 2 - x(l)x(5)x(6), x(l)x(6) 2 - x(0) 2 x(5), x(l) 2 x(6) - x(0) 3 , 

x(3) 2 x(4) - x(0)x(l) 2 , x(4)x(5) 2 - x(0)x(l)x(5), x(0)x(l)x(6) - x(2) 2 x(3), 

x(0)x(4)a;(5) - x(0) 2 x(l),x(0)x(5)x(6) -x(3)x(4) 2 , 

x(l)x(2)a;(6) - x(0)x(4)x(5), x(2) 4 - x(0) 4 , x(0) 3 x(5) - x(3) 3 x(4), 

x(0)x(5) 3 - x(3)x(4) 3 ,x(2) 3 x(3) - x(0) 3 x(l),x(2)x(3) 3 - x(0)x(l) 3 , 

x(0) 2 x(5) 2 - x(2) 2 x(3) 2 ,x(0) 2 x(l)x(5) - x(3) 2 x(4) 2 ,x(2) 2 x(3) 2 - x(0) 2 x(l) 2 , 

x(l) 2 x(2) 3 - x(0) 5 ,x(l) 4 x(2) 2 - x(0) 6 ,x(l) 6 x(2) - x(0) 7 ,x(l) 8 - x(0) 8 . 

This computation can be performed by applying the algorithm SigmaGBasis to 
the E-ideal L C P in the same way as for the example in Section 5. For details 
about different strategies to implement this method we refer to [18]. Note that 
T- ideals are also called "symmetric ideals" in [23] . 

9. Conclusions and further directions 

In this paper we show that a viable theory of Grobner bases exists for the alge- 
bra of partial difference polynomials which implies that one has symbolic (formal) 
computation for systems of partial difference equations. In fact, we prove that such 
Grobner bases can be computed in a finite number of steps when truncated with re- 
spect to an appropriate grading or when they contain elements with suitable linear 
leading monomials. Precisely, since the algebras of difference polynomials are free 
objects in the category of E-algebras where S is a monoid isomorphic to N r , we ob- 
tain the latter result as a Noetherianity criterion for finitely generated E-algebras. 
Among such Noetherian E-algebras one finds polynomial algebras in a finite num- 
ber of variables where a tensor product of a finite number of algebras generated 



NOETHERIAN QUOTIENTS OF THE ALGEBRA OF PARTIAL DIFFERENCE ... 19 

by single matrices acts over the subspace of linear forms. Considering that such 
commutative tensor algebras include group algebras of finite abelian groups one 
obtains that there exists a consistent Grobner basis theory for ideals of finitely gen- 
erated polynomial algebras that are invariant under such groups. In our opinion, 
this represents an interesting step in the direction of development of computational 
methods for ideals or algebras that are subject to group or algebra symmetries. 

As for further developments, we may suggest that the study of important struc- 
tures related to Grobner bases like Hilbert series and free resolutions should be 
developed in the perspective that their definition and computation has to be con- 
sistent to the symmetry one finds eventually on a polynomial algebra. An impor- 
tant work in this direction is contained in |15j . Finally, the problem of studying 
conditions providing S-Nocthcrianity (instead of simple Noetherianity) for finitely 
generated S-algebras is also an intriguing subject. 
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